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ABSTRACT 



A CDMA radiocommunication signals receiver for receiving l 
signals obtained from spectrum symbols spread using J 
pseudo -random sequences and having been propagated : 
along a number of paths. The receiver includes a filter 
configured to restore L unspread signals for each symbol, » 
corresponding to L different paths, a calculating circuit 
configured to calculate L estimates of the L different paths, I 
and a demodulator configured to process each of the L j\ 
unspread signals using the corresponding L estimates to » 
obtain L path contributions. Also included is an adder 
configured to form a sum of the Lpath contributions and for 
outputting an estimate of a received symbol, and a decision 
circuit configured to make a decision about a value of the 
received symbol based on a value of the estimate of the 
received symbol output by the adder. Further, the receiver 
processes blocks of N symbols, each block having data 
symbols and control symbols, each symbol being identified 
by a rank k that it occupies in the block, where k varies from 

0 to N-l. Also, for each path identified by an index 1, where 

1 varies from 0 to L-l, and for each block, the receiver 
considers a vector C, with N components that characterizes 
the path during the block, and the receiver defines a vector 
base B^, vectors of the vector base B^ being N eigenvectors 
of the matrix E [C/C/ 7 }, each vector C ; being decomposed 
in the vector base, where decomposition coefficients denoted 
Gut form independent random Gaussian variables. In 
addition, coefficients G^, define a vector G ; with N com- 
ponents for each path 1, and the calculating circuit estimates 
each vector G„ using an iterative process based on EM 
estimation-maximization algorithm based on a maximum a 
posteriori probability criterion. 

6 Claims, 5 Drawing Sheets 
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ITERATIVE RAKE RECEIVER AND 
CORRESPONDING RECEPTION PROCESS 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The purpose of this invention is an iterative rake receiver 
and a corresponding reception process. It is used in radio- 
communications applications and more particularly in the 
Code Division Multiple Access (CDMA) technique. It may 
be applied to the system defined by the IS' 95 standard and 
in third generation UMTS and IMT-2000 systems. 

2. Discussion of the Background 

In the CDMA technique, information symbols that are to 
be transmitted do not directly modulate a carrier, but are 
firstly multiplied by pseudo-random sequences (or codes) 
which have the effect of spreading their spectrum. On 
reception, the received signal is unspread by filtering 
adapted to the sequence used in transmission (or by 
correlation) and is then demodulated. 

This technique enables several users to use the same 
radiocommunications channel, provided that a particular 
sequence is assigned to each. 

Amultipath channel is generally used, in the sense that the 
radio-electric wave propagates along several different paths 
between the place of transmission and the place of reception. 
Therefore, the receiver does not receive a single signal for 
each transmitted information symbol, but it receives several 
more or less delayed and more or less modified copies. In 
order to reliably restore the transmitted information, the 
largest possible number of these copies has to be taken into 
account and they have to be recomposed in the receiver. 

One way of doing this was to design a special receiver, 
called a rake receiver, in the sense that it "rakes" information 
at different instants symbolizing the teeth of the rake. This 
type of receiver separates signals corresponding to different 
paths and includes several "teeth" or "branches" or "pins", 
each of which processes one of these signals. The signal is 
unspread and demodulated in each tooth. The next step is to 
recompose all signals in an adder. 

A rake receiver was initially described by R. PRICE and 
P, E. GREEN in an article entitled "A Communication 
Technique for Multipath Channels" published in the "Pro- 
ceedings of the IRE", vol 46, March 1958, pp 555-570. 

TTiere is also a description of this receiver in the general 
book by J. G. PROAKIS entitled "Digital 
Communications", 3 rd edition, MCGRAW-HILL, 1995, 
(third edition) 1989, (second edition). 

FIG. 1 attached diagrammaticaUy shows a receiver of this 
type. As shown, this receiver comprises a main input E, a 
filter 10 with a width adapted to the signal spreading band, 
L means 12 0 , . . . , 12/, . . . , 12^.! of restoring L signals 
unspread in frequency corresponding to L paths (these 
means usually include a filter adapted to one of the pseudo- 
random sequences used in transmission or a correlator and 
means of searching for signal peaks), L means 14 0 , . . . , 
14 ; , . . . , 14^_ 3 of estimating the characteristics of the L paths 
used by the various signals, L demodulation means 16 0 , . . . , 
16 /r . . . , 16 z-1 combining the unspread signals and estimate 
of paths, an adder 18 adding the N contributions output by 
the N demodulators, and finally a decision circuit 20 that 
outputs transmitted or reference symbols on a general output 
S in order to test the communication. 

After a strong demand from services requiring ever 
increasing throughputs, the frequency spreading band of 
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CDMA systems is continuously being increased. This 
increase in the band is accompanied by a continuous 
increase in the number of paths received at the receiver. For 
a given received power, this increase in the number of paths 

5 tends to reduce the power received per path and therefore to 
reduce the global quality of the channel estimate. 
Consequently, a constructive combination of contributions 
from these paths to the receiver is rarely guaranteed and can 
result in a severe loss of transmission quality. 

10 Since CDMA systems are inherently limited by multiple 
access interference, this loss of performance cannot be 
compensated by increasing the power. Furthermore, a solu- 
tion involving an increase in the number of pilot symbols 
would reduce the system capacity. 

15 The purpose of this invention is to correct this disadvan- 
tage. 

SUMMARY OF THE INVENTION 

The main purpose of this invention is to increase the 

2Q performances of CDMA systems by improving the reception 
quality for a given transmitted power, and therefore with an 
unchanged level of multiple access interference. This quality 
improvement can increase the capacity and coverage of the 
CDMA system. This improvement is obtained by optimizing 

25 operation of the receiver in the classical case of slow fading, 
but also in the more difficult case of very fast fading. 

Another purpose of the invention is to facilitate the 
construction of terminals by making them much less sensi- 
tive to inaccuracies in the local oscillator used to transpose 

30 the received signal in the base band. 

CDMA systems introduce the concept of a power control 
period (PCP). The power of the signal emitted by the 
transmitter remains constant during each of these periods but 
it can vary from one power control period to the next in order 

3S to compensate for slow fading (due to the distance and mask 
effects), and fast fading due to multipath effects when the 
terminal is moving slowly. Furthermore, control symbols are 
used in addition to symbols transporting information. The 
invention can reduce the relative number and/or power of 

40 these control symbols, for an equal reception quality. This 
objective is achieved by taking optimum account of control 
symbols for an arbitrary number of consecutive power 
control periods in the channel estimate. It is also achieved by 
taking account of control symbols (if any) included in the 

45 pilot channels of some CDMA systems in the optimal 
channel estimate, in the absence of matching antennas. It is 
also achieved by taking account of control symbols (if any) 
assigned to other users in the down link in the estimate, still 
in the absence of matching antennas. Finally, it is achieved 

S Q by taking optimum account in the estimate of this channel of 
all or some of the data symbols for these power control 
periods which are obviously more numerous and often 
contain more energy than the control symbols. The reduction 
in the number and/or the power of the control symbols is a 

55 way of consolidating the encoding of useful data and 
increasing the proportion of the transmitted power allocated 
to these data. 

With the invention, the optimum channel estimate can 
also take account of control symbols multiplexed in time and 
60 control symbols multiplexed on in-phase and/or 
in-quadrature components of the modulated signal. 

All power control period symbols can be used with the 
invention. Thus, shifts in the local oscillator frequency can 
be monitored and corrected, even if control symbols are 
65 grouped together. 

According to the invention, block by block processing is 
carried out every time that the received signal corresponding 
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to a given number of power control periods is available, like 
a conventional rake receiver, it always begins by unspread- 
ing the signals corresponding to significant paths selected 
for the final combination. It then makes an approximate 
estimate of the multipath channel by using control symbols 
associated with the received block only. This estimate is a 
way of characterizing the variation in the phase and the 
amplitude of each path during the block to be processed, for 
each symbol in the block, symbol by symbol. The receiver 
according to the invention demodulates and then combines 
the estimated path contributions and outputs a sample (or 
weighted output) for each data symbol contained in the 
block. 

In the case of a conventional rake receiver, these weighted 
outputs are used directly to detect and decode transmitted 
data symbols. These outputs have a certain reliability with 
respect to values taken on by the data symbols sent during 
a block. In the case of the receiver according to the 
invention, they may be used in addition to the control 
symbols, to provide an improved estimate of. each received 
path. This improved estimate of the multipath channel may 
be optimized, possibly tanking account of the encoded 
structure of the data symbols. Taking account of the cor- 
recting encoding results in better quality of the weighted 
outputs at the receiver. 

The weighted outputs obtained at the end of a given 
iteration may be used again, together with control symbols, 
to provide an additional improvement to the channel esti- 
mate. This improved estimate then improves the quality of 
the weighted outputs generated by the receiver. Therefore, 
the receiver output is looped back onto the estimating 
means. 

TTie optimum nature of the receiver according to the 
invention is related to the nature of the estimate of the 
multipath channel. This optimum nature depends firstly on 
the use of an iterative Estimation-Maximization (EM for 
short) type algorithm to construct the most likely channel, as 
a function of the received block. For example, this algorithm 
is described in the article by A. P. DEMPSTER, N. M 
LAIRD and D. B RUBIN entitled "Maximum Likelihood 
from Incomplete Data via the EM Algorithm", published in 
the Roy.Stat.Soc journal, Ser. 39, 1977. 

The optimum estimate for the channel also depends on the 
decomposition of each path received according to an expan- 
sion algorithm called the KARHUNEN-LOEVE algorithm. 
This decomposition enables flexible characterization of time 
variations of paths due to the oppler effect and can easily be 
included in the EM lgorithm itself. For example, the 
KARHUNEN-LOEVE lgorithm is described in the book by 
J. G. PROAKIS entioned above, 1989 version, pages 
340-344. 

' More precisely, the purpose of this invention is a CDMA 
radiocommunication signals receiver, these signals having 
been obtained from spectrum symbols spread using pseudo- 

» random sequences, these signal then having been propagated 
along a number of paths, this receiver comprising: 

means of restoring L unspread signals for each symbol, 

corresponding to L different paths, 
means of calculating L estimates of the L paths, 
demodulation means for processing each of the L 
unspread signals using the L corresponding estimates to 
obtain L path contributions, 
an adder to form the sum of these L contributions and to 

output an estimate of the received symbol, 
a decision circuit about the received symbol starting from 
the value of the estimate output by the adder, 



this receiver being characterized in that: 

a) it processes blocks of N symbols, each block 
comprising data symbols and control symbols, 
each symbol being identified by the rank k that it 

5 occupies in the block, where k varies from 0 to 

N-l, 

b) for each path identified by an index 1, where 1 
varies from 0 to L-l, and for each block, the 
receiver considers a vector C, with N components 

10 that characterizes the path during this block, 

c) the receiver comprises means of defining a vector 
base B^ these vectors being N eigenvectors of the 
matrix EfQC/ 7 }, each vector C, being decom- 
posed in this base, the decomposition coefficients 

15 denoted G lk forming independent random Gauss- 

ian variables, 

d) coefficients G lk , defining a vector G, with N 
components for each path 1, the estimating means 
being capable of estimating each vector G,, using 

20 an iterative process based on EM estimation- 

maximization algorithm based on a maximum a 
posteriori probability criterion. 
In a particular embodiment, the output from the adder is 
looped back onto the estimating means, the estimating 
25 means being used initially making use of control symbols 
contained in the block and assumed to be known, so that a 
first estimate of the data symbols contained in the block can 
be obtained from the output of the adder, the estimating 
means then using all estimated symbols present at the output 
30 from the adder and so on, the estimating means finally 
outputting the optimum value of G 7 (1=0, 1, . . . , L-l), after 
a final iteration. 

Another purpose of this invention is a process for recep- 
tion of CDMA type of radiocommunication signals, these 
35 signals having been obtained from spectrum symbols spread 
using pseudo -random sequences, these signals then having 
been propagated along a number of paths, this reception 
process comprising the following operations: 

for each symbol, L unspread signals corresponding to L 
40 different paths are restored, 

L estimates of the L paths are calculated, 
each of the L unspread signals is demodulated using the 
L corresponding estimates to obtain the L contributions 
45 of the paths, 

the sum of these L contributions is formed, which gives an 

estimate of the received symbol, 
a decision is made about the received symbol, based on 
the value of the estimate obtained, 
50 this process being characterized in that: 

a) blocks of N symbols are processed, each block c 
comprising data symbols and control symbols, 
each symbol being identified by rank k that it 
occupies within the block, where k varies from 0 

55 to N-l, 

b) for each path, identified by an index 1 where 1 
varies from 0 to L-l, and for each block, a vector 
C; with N components is considered, which char- 
acterizes the path during this block, 

60 c) the matrix E [QC/ 7 ] which has N eigenvectors 

denoted B* is considered, and these N eigenvec- 
tors B* are used as a base, each vector C, in this 
base is decomposed, the decomposition coeffi- 
cients denoted as G tk forming independent Gaus- 

65 sian random variables, 

d) the coefficients G /Jt defining a vector G, with N 
components for each path 1, and each vector G, is 
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estimated using an iterative process based on an 
estimation-maximization (EM) algorithm based 
on a maximum a posteriori probability criterion. 
In one particular embodiment, the iterative process is used 
initially taking account of the control symbols contained in 
the block and assumed to be known, which gives a first 
estimate for data symbols contained in the block, the said 
iterative process then taking account of all symbols in the 
block according to this first estimate, which can give a 
second estimate of the symbols in the block, and so on, until 
a satisfactory estimate is obtained for values of G,, that are 
used for the demodulation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1, already described, shows a known rake receiver; 

FIG. 2 illustrates the general structure of an iterative rake 
receiver according to the invention; 

FIG. 3 shows the estimating and optimization block 
according to the maximum a posteriori probability criterion; 

FIG. 4 compares the performances of a receiver according 
to the invention with two conventional rake receivers, giving 
the binary error rate as a function of the E/I 0 ratio, where I 0 
is the spectral power due to thermal noise and interference 
due to multiple accesses, and E is the average energy per 
received symbol; 

FIG. 5 compares the performances of the receiver accord- 
ing to the invention with two conventional receivers, giving 
the binary error rate as a function of the position of data 
symbols in a power control period for E/I 0 equal to 10 dB. 

DESCRIPTION OF PARTICULAR 
EMBODIMENTS 

The following mathematical developments, like the 
figures, are applicable to the case in which symbols are 
multiplexed in time only, for simplification and notation 
reasons. The case in which the symbols are multiplexed in 
time and by code on components in quadrature is dealt with 
in a similar manner. 

The receiver according to the invention searches for a 
representation of the multipath channel according to a 
criterion called the maximum a posteriori (MAP) probability 
criterion. This estimate requires knowledge of the emitted 
symbols, or at least their probability. In one particular 
embodiment, the receiver may use estimates present at the 
adder output. Therefore in this case, the output is looped 
back onto the various estimators as shown in FIG. 2. This 

I figure shows that the output from the adder 18, that outputs 
a signal denoted (where k denotes the rank of the symbol 
being processed, this rank varying from 0 to N-l as 
explained later) is looped back on estimators 14 0 , . . . 
14/, . . . , 14 z _j. But this loop back does not appear explicitly, 
as will become clear later. 

v The properties of the multipath channel are determined 
using a posteriori probabilities of some magnitudes, consid- 
ering the signal R lk (t) output by each unspreading circuit. 
This calculation is made by decomposing the signal into 
components, ensuring that the components are not corre- 
lated. This is done using the KARHUNEN-LOEVE decom- 
position algorithm. The principle of this algorithm is briefly 
described below. 

A random function of time z(t), with an average value 
equal to zero, has a correlation function denoted <KVc) equal 
to 
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5 where E is the mathematical expectancy. The function z(t) 
may be developed in series in the form: 

10 -i 

where the values of z„ are development coefficients and the 
values of f^,(t) are orthonormal functions over a given 
interval 0, T. 

Each coefficient z„ may be obtained from the function z(t) 
and the eigenfunctions f„(t) using the relation: 

20 It is demonstrated that the orthonormal functions f„(t) are 
the eigenfunctions of the equation: 

where the K m values are eigenvalues. 
25 If the signal output by an unspreading circuit (matched 
filter or correlator) is denoted R(t), we have: 

where C(t) is a function that characterizes the path taken by 
30 the wave, A(t) is the emitted information and N(t) is additive 
Gaussian noise. In order to estimate each path, the function 
C(t) that characterizes it will be decomposed using the 
KARHUNEN-LOEVE algorithm. In reality, CDMA type 
radiocommunications use digital magnitudes, or samples, 
35 rather than continuous functions in time. The output signal 
from the unspreading circuit for tooth order 1 in the rake is 
in the form: 

40 where k is the rank of the symbol, assumed to vary from 0 
to N-l if a block of N symbols is processed. For a given 
path, identified by the index 1, the N components C a form 
the N components of a vector C,. This is the vector that will 
be decomposed using the KARHUNEN-LOEVE algorithm 

45 which in this case involves discrete summations instead of 
integrals. Therefore, we will talk about "eigenvectors" rather 
than "eigenfunctions", but the spirit of the decomposition 
remains the same. 

Having mentioned these general aspects, the receiver and 

50 the process according to the invention will now be described 
in detail as follows, for the case in which a single power 
control period is considered. 

Each power level, corresponding to one power control 
period, is composed of N^ data symbols and N c control 

55 symbols, all modulated in phase (MDP2, MDP4, 
MDP8, . . . ). In the following, the total number of these 
symbols will be denoted N (N^N^+Nc), data symbols will 
be denoted ao, a lt . . . , ^ ND . 1 control symbols will be 
denoted a^ a^j, . . . , a^. 

60 These two categories of symbols may be spread using 
different spreading factors by means of finite pseudo- 
random sequences. Furthermore, they may be multiplexed 
either separately on in-phase and in -quadrature components 
of the modulated signal, or together in time. In the 

65 following, the relative position of each of these symbols in 
time with respect to the beginning of the corresponding 
power control period, will be. denoted P,-. 
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Id general, control symbols are composed firstly of N p 
pilot symbols a^, a^j, . . . , * ND _ NF -i known to the 
receiver and N C -N F symbols z„ D+NI » ^ ND ^ NP ^ - - • , 
dedicated to power control on the reciprocal channel and the 
indication of the received data format. 

During one given power control period, the emitted 
energy allocated to symbol a f will be denoted E,-. This energy 
is normally common to symbols in the same category, but it 
may differ in different categories. In this case, the energies 
assigned to the data and control symbols will be denoted Eq 
and E c respectively. 

The multipath channel seen by the emitted CDMA signals 
is composed of several paths which have or which could 
have temporal variations due to the Doppler effect. Each 
path is characterized by a given average power and a given 
Doppler power spectrum (DPS) that depend on the environ- 
ment and the speed of the mobile. Furthermore, fading on 
each path may be of the Rayleigh type or the Rice type. 

In general, Doppler power spectra with Rayleigh fading 
are either of the conventional type or the flat type. Conven- 
tional Doppler power spectra are encountered particularly in 
environments outside buildings, while flat Doppler power 
spectra are encountered particularly in environments inside 
buildings. 

The Doppler channel spreading effect is denoted B^,, and 
the Bessel function of the first type of order 0 is denoted 
J 0 (.). The self -correlation function of an average power path 
<|>(0) is then given by: 

in the case of a conventional DPS, and by 

<K*H>(0) sin (Aff D iT)/jifl 0 T, 

in the case of a flat DPS (Doppler power spectrum) The 
average power $(0) varies from one path to another and 
therefore characterizes the intensity profile of paths. 

As mentioned above, a rake receiver is composed of L 
teeth that can be used to follow the L paths with the highest 
power, and to constructively combine the contributions of 
these L paths. The number of paths used by a rake receiver 
is usually less than the real number of paths actually 
received. This number depends on the environment (external 
or internal) and the CDMA spreading factor. Typical values 
of 2 to 3 in internal environments and 4 to 8 in external 
environments are often used. 

Let R /Jt be the correlator output signal from the X th tooth 
corresponding to the K 1 * symbol a A transmitted during an 
arbitrary power control period. This signal may be in the 
form: 

where C lk is the gain factor of the X th path of the receiver seen 
by the symbol a^ and N tt is a complex noise that includes 
the thermal noise and the interference caused by multiple 
access of other mobiles. In order to simplify the analysis and 
design of the receiver, this noise is assumed to be Gaussian 
and decorrelated, and its variance is denoted I 0 . 

The gain factors for a given path are also assumed to be 
independent from each other. The reason for this is that two 
signals reaching the receiver with different delays have a 
good chance of not using the same paths and not meeting the 
same obstacles. Nevertheless, gain factors on a given path 
are usually correlated with each other. If E[.] and <(>,(.) denote 
the mathematical expectancy operator and the continuous 
self-correlation function of the 1 path respectively, then the 
discrete self -correlation function corresponding to this path 
is given by: 



4,740 Bl 

8 

where P ; and P- values identify the position in time of the 
symbols with respect to the flow of the power control period. 

5 In general, the receiver has a vague idea of the value of 
the Doppler spread B^, and of the shape of the Doppler 
power spectrum. Consequently, it adopts the representation 
of the least predicable multipath channel with a flat Doppler 
power spectrum, for which the Doppler spread is greater 

10 than or equal to the real spread. For simplification reasons, 
this upper limit of the spread will also be denoted B^,. It may 
be permanently fixed at the receiver as a function of its 
maximum authorized or actual speed. It may also be esti- 
mated by trial error, for example using symbols and/or the 

15 pilot channel. 

During each power control period, the receiver needs the 
most precise possible estimate of gain factors C Jk corre- 
sponding to data symbols, power control symbols and the 
format of the data received. Consequently, the receiver 

20 according to the invention is capable of taking account of the 
time correlation of gain factors of all paths combined by the 
receiver. It is also capable of taking account of all or part of 
the encoded structure of data and control symbols unknown 
to the receiver in order to improve this estimate. Finally, it 

25 can take account of all or some nearby power control period 
data and control symbols in order to optimize its estimate for 
a given power control period. 

For simplification purposes, we will only describe the 
case of a multipath channel estimator using symbols 

30 received during a power control period exclusively, in order 
to estimate the corresponding channel output. The transpo- 
sition operator will be denoted (.) r and the vector will be 
introduced: 

35 *HWW.J^-i) r 

The components of this vector are the N samples received 
during the power control period to be processed and corre- 
sponding to the 1* path selected by the receiver. The estimate 
40 of the multipath channel during this power control period is 
based entirely on these L vectors Rq, R 1} . . . , R^ of 
samples received. 

The modulus operator is denoted |.|. Remember that the 
amplitude |a J of a* that is equal to vl^ depends not only on 
45 the power control period, but also on the index of the emitted 
symbol. In practice, this amplitude is the same for each 
symbol category. 

The normalized transmitted symbols vector 

50 A-tAoAt Any- 

where A^-a^aJ, is introduced in order to eliminate this 
dependence on the amplitude. The k** component of the l** 
vector of received samples can therefore be in the form: 

55 

Rj^^CfjA t+Ntjk 

where C Ik is the k rt component of the 

60 C r Qa 0 \CMC n , . . . ,K_ 1 \C tJ ,_ 1 f 

vector of normalized gain factors corresponding to the I th 
path. This is the vector that is to be represented. 

This representation is based on the KARHUNEN- 
65 LOEVE decomposition algorithm, and consists of express- 
ing each of the L normalized vectors C 7 , 1=0, 1, . . . , L-l in 
the form: 
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Given the received vectors R„ the EM algorithm begins 

^ by calculating initial conditions G, (0) for the vectors G, 

Q = 2j °tk B k . starting from the received samples corresponding to the pilot 

symbols. 

5 The set of possible values taken on by the k** symbol of 

where the B k values are the N normalized eigenvectors of the a power period is denoted as s ^ the M of indexes 

covariance matrix F^EfCjC/ 7 ] of C, and the coefficients G tt of ilot symbols of a power ^niicl period are denoted s> and 

are independent random Gaussian variables with zero aver- ^ ^ taken 0fl b ^ ^ bo , A ^ aQ mdex 

age an ^ variance equal to the eigenvalues of the matrix F,, ^ s {& denoted ag M ^ of ^ 

remembering that the index k vanes from 0 to N-l. 10 . 4 . „ w . , iL _ * . . ., 

w. • i .1 , ,i i f4 , ^ , . initialization of the EM algorithm, the receiver has no idea 

It is assumed that the shapes of the Doppler power spectra , , , f j 4 « , « » ~ 

of the paths are identical and therefore that the correspond- abc f the val »f s °* data u s ^ ho]s A and the ^f°* T? / 

ing eigenvectors B* are also identical. If the exact charac- ™ to ™ ™nditional probability density P^R,^; 1 , 

teristics of the multipath channel are known, the (i j>* entry { G '~ " > for A * values that are not P Jot svmbols - If the 

in the matrix F will be: 15 modulation constellation has central symmetry (MPD2, 

7 MPD4, MPD8, . . . ) then the n" 1 component of the initial 

fuf^i^rPi)^^ condition of the X th path G/ 0) may be chosen to be equal to 

In practice, neither the shape of the Doppler power , 0 , _ , 

spectrum and the corresponding Doppler spread, nor the 2 o ^j^i * W * 

powers at which the power control periods are emitted, are 
known precisely. The receiver can then adopt the flat DPS 

with width B D representing an upper limit to the real where W, M is a weightingiactor equal to 

Doppler spread. Furthermore, it can assume that the average 

received power does not vary very much from one power 2 s 1 

control period to the next. In this more realistic case, the (i, W}a = i + / 0 y Cb 

)) th entry in the matrix F ? , given by 

fflrM 0 )^ Bin (xB D (P r fy)frB D (PrPj) This factor depends on both T^, the n th eigenvalue in the 

3 0 matrix F (which depends on the average power ty, (0) of the 

Do teTs 2d ° 0n UPPef '* Path and thC D ° Ppler Splead Bd ' Md ak ° th6 Cmitt6d 

1 F . " /. 4 . . « , - . energy associated with data and control symbols E„ and E c ) 

Also in practice, the receiver may have a bank of eigen- , ~r . . . , ,. At _ , . c < 

f r u * • i i c .i. i* % /\u and the noise vanance I 0 including the thermal noise and 

vectors for each typical value of the upper limit of the interference b multi le a 

Doppler spread in order to be able to better adapt to the m r erence v mu P e access - 

terminal speed. 35 Still based on all received vectors R,, the EM algorithm 

The channel estimator according to the invention makes then results in an iterative calculation of the reestimate G/* +1 

an iterative estimate of the multipath channel with fading, starting from the estimate Gf making use of the expression: 
using the MAP criterion. In its estimate, it may take account 
of the characteristics of the appropriate channel 
representation, and also the values of the pilot symbols and 40 
the encoded structure of unknown symbols (including data ~ In 

symbols). 

The MAP estimate of G,, where 1 represents all values 
from 0 to N-l, of an implementation of G, of the multipath 

channel with fading, is the value 45 defining the relation, between the n* component of the 



r -I 



(d+l)'' 1 reestimate of the X th vector of the appropriate repre- 

{Glw/^-arg max P{{G] l L ~ 1 y {R}t~ x sentation of channel G, (d) . 

which maximizes the a posteriori probability density ^ iterative estimate according to the invention of the 
P({G/}^ _1 , {R/}/_o i_1 )- According to the invention and 50 adequate representation G, may be carried out^a limited 
using the EM algorithm, a solution as close as is desired to number D times such that the global estimate G/ 0) obtained 
the exact solution can be obtained iteratively. guarantees an imperceptible degradation of receiver perfor- 
in the general case, the a posteriori probability density mances compared with the optimum solution G,. 
P({ G /}/-o _1 » { R /L-o i_1 ) t0 be maximized has several global These operations are illustrated in FIG. 3, which shows 
maxima that lead to an ambiguity in the estimate of the 55 mc estimator rank P, namely 14,, with means 30 defining the 
channel using the MAP criterion. This ambiguity may be vectors base B k and means 32 for calculating weighting 
eliminated by using pilot symbols known to the receiver. coefficients. The estimator 14„ uses a first estimate repre- 
However, this is often insufficient, since this probability sentcd symbolicaUy by block L that outputs the initial 
density has also local maxima that may be reached by the estimate G/°>, and then an estimate of order d represented by 
EM algorithm instead of the single g obal maximum. This 60 ^ { ^ ^ Q «, men m d+1 med 

problem can be solved by using pilot symbols to make a , ? T A .\ ' J£ « f. 

r , t . . , 4 . 4 . • i • • . . „(o) , n by block \.. that outputs G/ a+i } , and finally a final iteration 

suitable determination of the initial conditions G, w , 1=0. t , J?* 1 .•/,,, , ^ 

I l^l of order D represented by block \ D that outputs G/ \ 

1 The ' EM algorithm reestimates the G, vectors by Therefore, the receiver that has just been described corn- 
induction, to guarantee monotonous growth of the condi- 65 prises L circuits like the circuit shown as reference 14,, in 
tional a posteriori probability density p({G,},_ 0 i ~ 1 , FIG. 3. Therefore, in order to refine its estimate of the G, (rf5 
{R,},^ -1 ). channel, the circuit 14, must have probabilities 
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in other words the probability that the symbol A* will be 5 
equal to one among all the possible values, considering R, 
and G,. 

We have seen that D+ 1 iterations are carried out according 

to the invention in order to calculate G (0 \ G<*> . . . , G^> in where tann t-] * the hyperbolic tangent function and *e {.} 

sequence. However, a problem arises for the first pass (d-0) 10 is the "real part" function. The arguments of the 9U {.} 

since at this time there is no available value of G,. Therefore operator are naturally supplied by the receiver during each 

it is impossible to calculate the probabilities defined above iteration in the estimating algorithm. At the end of the 

in a rigorous manner. v iterative process (as we will see below), these arguments 

u ?K !° t Tr n{1 ° n> ^ Pll0t ( ° r r k eferenc *> will be supplied directly by the receiver to an encoder (if 

bols for which probabilities are known, are then used. For a 15 x . . Li , , ..... 

pilot symbol the probability that A, will be equal to the value USin * wei S hted decodlD S at the 1D P ut 

D k is equal to 1 (and the probability that A k is not equal to The complexity of the channel estimator may be further 

the value D* is equal to zero). Equally shared probabilities reduced if one of the following clipping functions is used: 

will be used for other symbols (for example Vi and X A will be 

used for the two possible values for a binary symbol). 20 f * if -i s * s l 

For subsequent iterations (d*0), values of G/* (and R 7 ) *o(*) = 

are available through knowledge of A^ values outputs 
from the recombination of all demodulators defined by the 
expression 



x elsewhere 



^\ (N-\ 



25 " x 

d d>t B t ) f 0 if -1/2 *xsin 

*i<*) = 

[sign of x elsewhere 

30 



and therefore the probabilities can be calculated as far as the 
(D+l)* iteration (d-D). 

This process is rigorously applicable to cases in which 

symbols A* are information symbols not protected by a <D 2 (*)«sign of x 

correction code, as shown below (or possibly protected by a « . , „ „ 

code but which we do not want to use for reasons of ***** of ^ hyperbolic tangent function, 

simplicity). On the other hand, if the A* values are symbols If all steps in the EM algorithm are to be executed 

protected by a correction code that we do want to use, an correctly, the receiver must have an upper limit B D for real 

algorithm called the BAHL algorithm is used to obtain the Doppler spread, and an estimate of the noise variance ^ and 

probabilities P(A*=A|{R,};_ 0 Z_ \ {G,},^" 1 necessary for 40 the individual power of each selected path O/0) 

iteration d+1 starting from iteration d. This algorithm carries , f the rece iver has a bank of eigenvectors for different 

out a complex operation, which is actually an improved vahies of the Doppler spfead uppef ^ B it cm 

demodulauon process. Therefore the demodulaUon opera- { ad te base yectors R as a q{ ^ 

turn carried out by demodulators 16 0 , . 16, . . . , 16, ,s ^ of ^ D kf ea &nd ^ yM 

no longer in the loop since another operation, more complex 45 Q ^ ^ e terminal 
and more complete is already included in the iteration. In 

this case there is no longer any need to loop back on If & Q u PP er limit of lhe Doppler spread is not updated as 

estimators at the output from adder 18. a function of the real speed of the terminal, a unique value 

BAHL's algorithm is described by L. R BAHL, J. w* 11 be assigned to it and the receiver only outputs the 

COCKE, F. JELINEK and J. RAV1V in the article entitled 50 corresponding base vectors. 

"Optimal Decoding of Linear Codes for Minimizing Symbol Similarly, due to the quantities mentioned above, the 

Error Rate", published in IEEE Transactions on Information receiver can calculate and output the necessary weighting. 

Theory, vol. IT-2°, March 1974. Having produced a D rt estimate 0, (o) , this final estimate 

In the special case of repetition encoding, the encoded fa used as a representation of lhe cbanncl and this represen . 

structure may immedia.ejy be integrated into the expl.cit 55 u1ioa is denotcd G For each M of rank k> thc 

formula givmg the (d+l)" iteration, thus avoiding the need , , , . 1jC c ' , , , t D , 4 , . 

for the BAHL al orithm demodulator 16 ? forms the product of K lk by the complex 

Inonespecialfmb^dtoent.thecomplexityofthechannel quantity of ^ namely C',„ and the adder 18 

estimator may be reduced without any significant loss of forms ** f m of 31 ~° ,n ^ , K 10,,S from L paths and 

performance, by only selecting the eigenvectors correspond- 60 o^P^ a final signal A,< 'defined by: 
ing to important eigenvalues, in the channel representation. 

In another particular embodiment, the expression G / „ (<, * 1) a <d> g k c 

may be further simplified in a special case of modulations * ^ * u 
MDP2, MDP4, etc. For modulation MDP2, the symbol A k 

with index k will be equal to the values in the set S* 65 

including the two values +VE^ and -VE^ and the general If the values of G* ik are expressed as a function of the 

expression for G /M (rf+1) is transformed into basic vectors, we get 
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The A k ^ signals may be used by a Viterbi detector/ 
decoder to retrieve unknown symbols (data) transmitted 
during a power control period, when encoding was applied 
to protect the data. 10 

For MDP2 demodulation, the real parts of the A k ^ 
signals are sufficient for decoding and act as weighted 
.outputs. 

Furthermore, for an uncoded MDP2 modulation with 
^unknown symbols taking on the values +VEJ[ and -VE^ with 15 
[lequal probabilities, the decision on the symbol A* is simply 
[Jgiven by 

i4 4 -eign of BefA^} 

The invention that has been described above may be 20 
implemented regardless of the distribution of reference 
symbols. These reference symbols may be grouped or dis- 
tributed in any manner. In particular, this invention can be 
used with a particular distribution as described and claimed 
in the French patent application deposited by this applicant 25 
on the same day as this patent application, and entitled 
"CDMA digital communication processes with distribution 
of reference symbols". 

Operation of the receiver according to the invention has 
been simulated in service for the 8 kb/s channel on the up 30 
link of the UMTS's CDMA system. Data and control 
symbols are multiplexed on components of the transmitted 
modulated signal, in phase (I) and in quadrature (Q) respec- 
tively. 

The control symbol spreading factor was twice the data 35 
symbol spreading factor. The duration of one power control 
period was 0.625 ms, and contained 1^=20 data symbols 
with period ^=31 .25 /is and N c =10 control symbols with 
period T c *62.5 ^s- The time positions of the control sym- 
bols are given by: 40 

T Dt M), 1 , . . . , No-1 

for data symbols, and by 

for control symbols. Control symbols included N^=6 pilot 
symbols assumed to be grouped at the beginning of each 
power control period. Furthermore, the average power of the 
data symbols was assumed to be twice the average power of 50 
control symbols. Therefore the average received energy E is 
identical for all received symbols and satisfies the equation 



£= (Z /,(0) ) £i * Jt = 0 ' 1 



The terminal is assumed to move at a speed of 500 km/h 
and uses a carrier frequency of 1.92 GHz. The Doppler 
spread corresponding to these choices is 1.778 kHz. The 60 
channel is assumed to have L-3 paths with equal average 
power. 

The performances of a receiver according to the invention 
and satisfying these assumptions were compared with the 
performances of two conventional receivers used in CDMA 65 
receivers. The first uses the estimating algorithm based on 
averaging, which compensates for the modulation applied to 



received samples corresponding to pilot symbols, and uses 
the average of these values to estimate the multipath chan- 
nel. The second uses the linear estimating algorithm that also 
compensates modulation of pilot symbol samples, but car- 
ries out linear interpolation and/or extrapolation of the 
channel using the minimum mean square error (MMSE) 
criterion. 

The comparison is made by the variation in the unproc- 
essed binary error rate (BER) (ignoring possible error cor- 
rection encoding) for distributed or grouped pilot symbols, 
using the three proposed channel estimating algorithms 
(according to the invention, by averaging or by linear 
estimate). 

FIG. 4 shows the binary error rate as a function of El 0 (the 
ratio between the average energy received per symbol and 
the noise level), and in FIG. 5, the same rate is shown as a 
function of the position of data symbols within a power 
control period for E/I 0 =10 dB. The references on the curves 
in these two figures correspond to the following character- 
istics: 

40, 50: estimate by averaging, grouped pilot symbols, 

41, 51: estimate by averaging, distributed pilot symbols, 

42, 52: linear estimate, grouped pilot symbols, 

43, 53: estimate according to the invention, grouped pilot 
symbols, 

44, 54: linear estimate, distributed pilot symbols, 

45, 55: estimate according to the invention, distributed 
pilot symbols, 

46: theoretical curve, perfectly known channel. 

These results show that the receiver according to the 
invention always provides better performances than either of 
the other two conventional receivers. The two curves closest 
to the theoretical curve correspond to the invention. 

For example, for a gross BER equal to 2.10" 2 and grouped 
pilot symbols, the receiver according to the invention guar- 
antees a gain in El 0 of the order of 3 dB compared with the 
best conventional receivers. 

What is claimed is: 

1. A CDMA radiocommunication signals receiver for 
receiving signals obtained from spectrum symbols spread 
using pseudo -random sequences and having been propa- 
gated along a number of paths, said receiver comprising: 
means for restoring L unspread signals for each symbol, 

corresponding to L different paths; 
means for calculating L estimates of the L different paths; 
demodulation means for processing each of the L 
unspread signals using the corresponding L estimates to 
obtain L path contributions; 
an adder for forming a sum of said L path contributions 

and for outputting an estimate of a received symbol; 
a decision circuit for making a decision about a value of 
the received symbol based on a value of the estimate of 
the received symbol output by the adder, wherein 
the receiver processes blocks of N symbols, each block 
comprising data symbols and control symbols, each 
symbol being identified by a rank k that it occupies 
in the block, where k varies from 0 to N-l, 
for each path identified by an index 1, where 1 varies 
from 0 to L-l, and for each block, the receiver 
considers a vector C, with N components that char- 
acterizes the path during the block, 
the receiver further comprises means for defining a 
vector base B^, vectors of the vector base B^ being 
N eigenvectors of the matrix E [CjC/ 7 }, each vector 
Q being decomposed in said vector base B K , where 
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decomposition coefficients denoted form inde- 
pendent random Gaussian variables, and 
wherein coefficients G^, define a vector G, with N 
components for each path 1, and the calculating 
means estimates each vector G,, using an iterative 5 
process based on EM estimation-maximization algo- 
rithm based on a maximum a posteriori probability 
criterion. 

2. The receiver according to claim 1, wherein 

an output from the adder is looped back onto the calcu- 10 
lating means, and 

the calculating means are initially used based on control 
symbols contained in the block and assumed to be 
known, which gives a first estimate for data symbols 
contained in the block at the output from the adder, and 15 
the calculating means then uses estimated symbols 
present at the output from the adder and outputs a 
quasi-optimal value of the vector G„ after a final 
iteration. 

3. A process for receiving CDMA radiocommunication 20 
signals obtained from spectrum symbols spread using 
pseudo-random sequences and having been propagated 
along a number of paths, said process comprising: 

restoring for each symbol, L unspread signals correspond- 25 
ing to L different paths; 

calculating L estimates of the L different paths; 

demodulating each of the L unspread signals using the 
corresponding L estimates to obtain L contributions of 
the paths; 30 

forming a sum of said L contributions, which gives an 
estimate of a received symbol; 

making a decision about a value of the received symbol 
based on a value of the estimate of the received symbol; 

processing blocks of N symbols, each block comprising 
data symbols and control symbols, each symbol being 
identified by a rank k that it occupies within the block, 
where k varies from 0 to N-l; 

considering for each path, identified by an index 1 where 40 
1 varies from 0 to L-l, and for each block, a vector Q 
with N components, which characterizes the path dur- 
ing the block; and 

considering a matrix E [QC/ 7 ] which has N eigenvectors 
denoted B k using said N eigenvectors B* as a base, and 45 
decomposing each vector C t in said base, where decom- 
position coefficients denoted G lk form independent 
Gaussian random variables, 

wherein the decomposition coefficients G tk define a vector 
G, with N components for each path, and each vector 50 
G, is estimated using an iterative process based on an 
estimation-maximization (EM) algorithm based on a 
maximum a posteriori probability criterion. 

4. The process according to claim 3, wherein the iterative 
process is initially used based on control symbols contained 55 
in the block and assumed to be known, which gives a first 
estimate for data symbols contained in the block, and the 
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iterative process then takes into account all symbols in the 
block according to said first estimate, which gives a second 
estimate of the symbols in the block, until a satisfactory 
estimate is obtained for values of G,, that are used for the 
demodulation process. 

5. A CDMA radiocommunication signals receiver for 
receiving signals obtained from spectrum symbols spread 
using pseudo -random sequences and having been propa- 
gated along a number of paths, said receiver comprising: 

a filter configured to restore L unspread signals for each 

symbol, corresponding to L different paths; 
a calculating circuit configured to calculate L estimates of 

the L different paths; 
a demodulator configured to process each of the L 
unspread signals using the corresponding L estimates to 
obtain L path contributions; 
an adder configured to form a sum of said L path contri- 
butions and for outputting an estimate of a received 
symbol; 

a decision circuit configured to make a decision about a 
value of the received symbol based on a value of the 
estimate of the received symbol output by the adder, 
wherein 

the receiver processes blocks of N symbols, each block 
comprising data symbols and control symbols, each 
symbol being identified by a rank k that it occupies 
in the block, where k varies from 0 to N-l, 
for each path identified by an index 1, where 1 varies 
from 0 to L-l, and for each block, the receiver 
considers a vector C, with N components that char- 
acterizes the path during the block, 
the receiver defines a vector base B^, vectors of the 
vector base being N eigenvectors of the matrix E 
[CjC/* 7 }, each vector C, being decomposed in said 
vector base, where decomposition coefficients 
denoted form independent random Gaussian 
variables, and 
wherein coefficients G^, define a vector G, with N 
components for each path 1, and the calculating 
circuit estimates each vector G„ using an iterative 
process based on EM estimation-maximization algo- 
rithm based on a maximum a posteriori probability 
criterion. 

6. The receiver according to claim 1, wherein 
an output from the adder is looped back onto the calcu- 
lating circuit, and 

the calculating circuit is initially used based on control 
symbols contained in the block and assumed to be 
known, which gives a first estimate for data symbols 
contained in the block at the output from the adder, and 
the calculating circuit then uses estimated symbols 
present at the output from the adder and outputs a 
quasi-optimal value of the vector G,, after a final 
iteration. 
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